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ABSTRACT 

We investigate the clustering of galaxy groups and clusters in the SDSS using the iBerlind et alJ 
( 2006) group sample, which is designed to identify galaxy systems that each occupy a single dark 
matter halo. We estimate group masses from their abundances, and measure their relative large-scale 
bias as a function of mass. Our measurements are in agreement with the theoretical halo bias function, 
given a standard ACDM cosmological model, and they tend to favor a low value of the power spectrum 
amplitude a$. We search for a residual dependence of clustering on other group properties at fixed 
mass, and find the strongest signal for central galaxy color in high mass groups. Massive groups with 
less red central galaxies are more biased on large scales than similar mass groups with redder central 
galaxies. We show that this effect is unlikely to be caused by errors in our mass estimates, and is 
most likely observational evidence of recent theoretical findings that halo bias depends on a "second 
parameter" other than mass, such as age or concentration. To compare with the data, we study the 
bias of massive halos in N-body simulations and quantify the strength of the relation between halo 
bias and concentration at fixed mass. In addition to confirming a non-trivial prediction of the ACDM 
cosmological model, these results have important implications for the role that environment plays in 
shaping galaxy properties. 

Subject headings: cosmology: large-scale structure of universe — galaxies: clusters 



1. INTRODUCTION 

Our current understanding of galaxy formation is cen- 
tered on the well-motivated assumption that all galaxies 
are formed and live their lives within dark matter halos. 1 
The observed clustering of galaxies is thus tied to the 
clustering of halos. Indeed, the modern way of modeling 
galaxy clustering is to combine halo profiles, abundances, 
and clustering (measured from cosmological N-body sim- 
ulations) with prescriptions (usually referred to as the 
halo occupation distrib ution) that specify how galax- 
ies o ccupy halos (e.g. ISeliakI 1200(1 IPeacock fc SmithI 
2000HScoccimarro et all200lUBerhnd fc Weinber J2002I: 
Coorav fc Shethl 12002'). These halo models usually as- 
sume that the large-scale clustering of halos depends only 
on halo mass. Under this assumption, galaxy properties 
only correlate with large-scale environment via their cor- 
relation with halo mass. 

The dependence of halo clustering on mass has 
been studied extensively in the context of ACDM cos- 
mological models and is fairly well understood (e.g., 
i Cole fe Kaiserll989UMo fc Whiteil996llSheth fc Tormenl 
ll999HSeliak fe Warrenll2004t iTinker et al.ll2005(l . More- 
over, the exact form of this dependence is sensitive to 
the values of cosmological parameters, such as the mass 
density f2 m and the amplitude of the matter power spec- 
trum (78, and can thus, in principle, be used to constrain 
cosmological models. 

Recent studies using N-body simulations have shown, 
however, that the clustering of halos also depends 
on a second parameter, in addition to halo mass. 

1 Throughout this paper, we use the term "halo" to refer to 
a gravitationally bound structure with overdensity p/p ~ 200, so 
an occupied halo may host a single luminous galaxy, a group of 
galaxies, or a cluster. Higher overdensity concentrations around 
individual galaxies within a group or cluster constitute, in this 
terminology, halo substructure, or "sub-halos". 



iSheth fc Tormenl l|2004fl first noted this when they de- 
tected that halos in dense environments formed at 
slightly earlier times than hal os of the same ma ss that are 
in low density envir onments. iGao et aT] (|2005T ) and then 
lHarker et all 12006') explored this in a larger simulation 
and found a much stronger signal: early forming halos 
cluster more strongly on large scales than late-forming 
halos of the same mass. This effect, sometimes referred to 
as "assembly bias" , was strongest for halos of mass much 
less than the n onlinear mass M*, and it disappeared for 
larger masses. iWechsler et"aTl (|2006fl did a more thor- 
ough anal ysis of these iss ues and, in addition to confirm- 
ing the IGao et alJ l)2005|) result at low masses, found the 
opposite effect for high mass halos, with a turnaround 
at M ~ M*. High mass halos that have recently assem- 
bled c luster more strongly t han old halos of the same 
mass. IWechsler et alJ {2006) also showed that all these 
effects persist when using halo concentration, rather than 
age, as the "second parameter". This is not surprising 
given that halo concentrat ion corre lates with formation 
time <|Wechsler et al.ll200l . Finally. IWetzel et all (2006) 
studied the high mass regime in detail and confirmed the 
IWechsler et all l|2006j) results. 

The correlation of large-scale environment with halo 
age/concentration at fixed mass will propagate into a 
correlation with galaxy properties if these properties 
thems elves correlate with h alo age/concentration at fixed 
mass. iCroton et alJ ()2006j) studied this issue in a semi- 
analytic model and found that the clustering of red 
and blue model galaxies on large scales cannot be ex- 
plained solely by the color-mass and mass-clustering re- 
lations. In other words, galaxy color must also cor- 
relate with large-scale environment through a "second 
parameter" like halo age or concentration. This may 
not be surprising, but observational analyses have shown 
that, in the real universe, if these effects exist they 
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are quite small iAbbas fc Shethl200llSkibba et al.l2006t 
iBlanton fc Berrindll2006lK The apparent discrepancy be- 
tween the theoretical and observational studies could be 
due to the fact that the predicted effects for halos are 
only strong for either (1) very low mass halos that do 
not contain galaxies luminous enough to be included in 
the observational samples, or (2) very high mass halos 
that are rare and do not contribute much to the total ob- 
served clustering of galaxies. Alternatively, it is possible 
that the trends seen for halos do n ot, in fact, propagate 
into trends for galaxies and that the lCroton et al.| l)2006f) 
results do not apply to the real universe. It is worth not- 
ing that the semi-analytic model used in this study does 
not predict the correc t clustering of red and blue galaxies 
iSpringel et alJl200l . 

A more direct way of observationally probing the the- 
oretical results is to study not galaxies, but groups and 
clusters of galaxies, because these are the objects that 
presumably correspond one-to-one with halos. Given 
a galaxy luminosity limit near L„, which is typical for 
current large surveys, low mass halos (M ~ 1O 12 M0) 
will contain one isolated galaxy, intermediate mass ha- 
los (M ~ 1O 13 M ) will contain a group of galaxies, and 
high mass halos (M ~ 10 14 M Q ) will contain a cluster. 2 
Studying the clustering of groups can thus serve as a 
proxy for studying halos. In particular, we want to look 
for an observed group property that acts as a "second 
parameter" and shows a clustering dependence at fixed 
mass. Many studies have focused on measuring the clus- 
tering of groups as a function of properties that act as a 
proxy for mass ([Ba hc^Ilfc_Sj3neir^[1983; |GiuTicin_jrt_aL| 
| 200C||Bahcall et , alJl2003tlPadilla et alJl2004 lYang et all 
120051 iCoil et al.ll2006fl . but few have looked at other pa - 
rameters. One notable exception is lYang et alJ (|200(tT| . 
who found evidence of such an effect using a group cat- 
alog derived from the Two Degree Field Galaxy Red- 
shift Survey (2dFGRS; IColless et alJl2?X)T^ . Specifically, 
they found that, at fixed estimated group mass, groups 
containing central galaxies with low star formation rates 
(SFR) are more clustered on large scales t han groups 
contai ning central galaxies with high SFR. lYang et alJ 
( 2006) found this result for all group masses in their sam- 
ple. 

In this paper we investigate the clus tering of groups i n 
the Sloan Digi tal Sky Su r vey (S PSS. lYork et all 1200^ . 
We use the iBerlind et alJ l)2006l) group catalog that was 
constructed using an algorithm tuned to identify galaxy 
systems that each occupy a single DM halo. We describe 
the galaxy and groups samples that we use, as well as our 
analysis methods in § El We first study group clustering 
as a function of estimated mass in § and then search 
for a residual dependence on other group properties in 
§ EI In § we analyze the clustering of halos in N-body 
simulations to compare with our group measurements. 
We discuss and summarize our results in §Eland[3 

2. DATA AND ANALYSIS 
2.1. SDSS 

The SDSS is a large imaging and spectroscopic sur- 
vey that is mapping two-fifths of the Northern Galactic 
sky and a smaller area of the Southern Galactic sky, us- 

2 Throughout this paper, we use the term "groups" to describe 
the observational counterparts of halos in all three of these regimes. 



ing a dedicated 2.5 meter telescope l|Gunn "eTalll2006T) at 
Apache Poi nt, New Mexico. T he survey uses a photomet- 
ric camera ijGunn et al]ll99cf) to scan t he sky simulta- 
neously in five photome tric bandpasses ijFukugita et alJ 
119961 ISmith et alJl2002]) down to a limiting r-band mag- 
nitude of ~ 22.5. The imaging data are processed 
by au tomatic software that does astrometry ijPier et alJ 
2003), source identi f ication, deble nding and photometry 
lILupton et al.ll20011 ILuptonll2005lh p hotometric calibra- 
tion jHogg et all 120011 ISmith et all 120021 iTucker et alJ 
2006), and data quality assessment ijlvezic et all 12004 K 
Algorithms are applied to select spectroscopic targets for 
the main galaxy sample l|Strauss et al.ll2 002h the lumi- 
nous red galax y sample |E iseTisteT n^t^lTl2001|) . and the 
quasar sample l|Richards^t"1i!"ll2002[) . The main galaxy 
sample is approximately complete down to an apparent 
r-band Petrosian magnitude limit of < 17.77. Targets 
are assigned to spectroscopic plates using an adaptive 
tiling algorithm ijBlanton et al.H2003dj) . Finally, spectro- 
scopic data reduction pipelines produce galaxy spectra 
and redshifts. 

Our group sample is identified from the large-scale 
structure sample sample 14 from the NYU Value Adde d 
Galaxy Catalog (NYU-VAGC; IBlanton et all H^). 
Gala xy magnitudes are corrected for Galactic extinc- 
tion iSchlege l et alj|1 998[) and ab solute magnitudes are 
k-corrected IIBlan ton et al.||2003al) and corrected for pas- 
sive evolution ( Blanton et al. 2003c) to rest-frame magni- 
tudes at redshift z = 0.1. The galaxy sample that we use 
was made publicly av ailable (and superseded) with th e 
SDSS Data Release 4 i jAdelman-McCarthv et alll2006|) . 
We restrict our sample to regions of the sky where the 
completeness (ratio of obtained redshifts to spectroscopic 
targets) is greater than 90%. Our final sample covers 
3495.1 square degrees on the sky and contains 298729 
galaxies. 

2.2. Group and Cluster Samples 

We use the group sample described in IBerlind et alJ 
(2006). Groups were identifi ed using a redshi f t-spac e 
friends-of- friends algorithm (Gcllcr & Huchra 1983), 
which was applied to a volume-limited sample of galaxies 
spanning the redshift range from z — 0.015 to 0.1 and 
complete down to an absolute r-band mag nitude of -19.9 
(the Mr20 sample in IBerlind et al.l l2006K This sample 
contains 57138 galaxies. With the help of mock galaxy 
catalogs, the group-finding algorithm was tuned to iden- 
tify galaxy systems that each occupy the same underlying 
dark matter halo. The result was a group sample that 
has unbiased richness and size distributions with respect 
to the underlying halo population. We emphasize that 
the group-finding algorithm used only the galaxy posi- 
tions in redshift space and did not use other properties, 
such as color. The volume-limited sample, mock cata- 
logs, group-finding algorithm, a nd the resulting gro up 
catalog are described in detail in Bcrlind et alJ ((2006) . 

We compute a total r-band luminosity for each group 
by summing the luminosities of its member galaxies. 
These group luminosities are not actually "total", since 
they do not include the light of member galaxies below 
the -19.9 absolute magnitude threshold. However, since 
the groups were identified from a volume-limited sample 
of galaxies, these luminosities should roughly preserve 
the rank order of true group luminosities. We then ob- 
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tain masses for our groups by matching the group lu- 
minosity function to a theoretical halo mass function, 
assuming a monotonic relation bet ween group luminos- 
ity and mass. Specifically, we use a iWarren et all pOOl 
halo mass function with the following values for cosmo- 
logical parameters: fl m — 0.3, = 0.04, n s = 1.0, 
as = 0.9. We use this slightly outdated cosmological 
model because the mock catalogs that were used to test 
the group-finder were based on this model. We thus ob- 
tain group mass esti mates from their abundances, which 
iBerlind et alJ l)2006|) showed to be unbiased with respect 
to the abundances of halos. Our mass estimates ig- 
nore the scatter in the mass-luminosity relation, but this 
should not affect our results much because we only use 
the group masses to create a set of mass-threshold group 
samples. However, we note that these are actually group 
luminosity threshold samples, which are characterized as 
having the same abundance as halos of the stated mass 
thresholds. 

We create four group sub-samples with the following 
mass thresholds: log(M/ft- 1 M ) > 14.0, 13.5, 13.0, 
and 12.5. The resulting number of groups in each of 
these samples is 327, 1316, 4299, and 12655, respec- 
tively. Approximately 45% of groups in the mass range 
12.5 < log(M / 7i _1 Mq) < 13.0 are isolated galaxies (i.e., 
N = 1 groups). On the high mass end, ~ 90% of groups 
more massive than lag(M/h M®) > 14.0 contain 9 or 
more galaxies. 

We define the position of each group on the sky to 
be the centroid of its member galaxy positions, and the 
group redshift to be the mean member galaxy redshift. 
Our results are not sensitive to these choices, since we 
focus on the clustering of groups on large scales. 

2.3. Analysis Method 

We wish to measure the relative bias of groups as 
a function of mass and other group properties. The 
bias between two samples A and B is typically defined 
as the large-scale asymptotic value of ^£,AA(r)/£,BB(r), 
where £,aa{t) and £-Bs( r ) are the two-point autocorrela- 
tion functions of the two samples. However, our group 
samples are small enough (especially at high masses) that 
the autocorrelations are very noisy, since the number of 
group-group pairs is small. We therefore measure cross- 
correlations of our group samples with the full volume- 
limited galaxy sample from which the groups were iden- 
tified. These cross-correlation functions are much higher 
signal-to-noise because the number of galaxies is much 
higher than the number of groups. The bias between 
group samples A and B is now the large-scale asymp- 
totic value of b(r) = £ag(?) / '^bgM , where £,AG( r ) and 
£bg i r ) are the cross-correlation functions of group sam- 
ples A and B with the galaxy sample G. This bias will be 
identical to the autocorrelation bias, as long as the cross- 
correlation coefficients of the two group samples are the 
same. 

We estimate the cross- correlation functi o n usin g a sym- 
metrized version of the ILandv fc Szakwl JT993) estima- 
tor, 

RR \nin 2 J RR \ n\ J RR \n 2 ) 

(1) 



where B\Di are data-data pair counts between the 
two samples (in this case groups and galaxies), D\R and 
D 2 R are data-random pair counts between each sample 
and a catalog of random points, and RR are random- 
random pair counts within the random catalog, ni, n 2 , 
and tir are the number densities of group, galaxy, and 
random catalogs. Since the group and galaxy samples 
are both volume-limited and have the same geometry, it 
is sufficient to use one random catalog for both. We use 
a random catalog containing one million points, which is 
large enough that the noise in the estimated correlation 
function is dominated by the data-data term. 

Since our data are in rcdshift-space, where galaxy and 
group peculiar velocities distort their positions along the 
line-of-sight, we bin all pairs in a grid of perpendicular 
and line-of-sight separations: r„ and it. We est imate 
these separations according to iFisher et al.l (|l994) . For 
two points with redshift positions vi and V2, we define 
the vectors s = vi — V2 and 1 = g(vi + V2). r p and it for 

these points are then tt = (s ■ 1)/|1| and r p = %/s • s — n 2 . 
Once we estimate Ci2(?>, tt), we integrate along the red- 
shift direction to get t he projected correlation function 
Ipavis & Peebleslll983l) . 

■w P ,i2(r P ) = 2 dn £ 12 (r p ,ir). (2) 
Jo 

In practi ce, we only in t egrate out to 7r max = 40/i _1 Mpc, 
following IZehavi et all (|2002j) . 

Once we have measured the projected group-galaxy 
cross-correlation functions for two group samples A and 
B, we define a bias function 

b(r p ) = w p ,AG(rp)/w Pt BG(r P )- (3) 

We then compute the asymptotic large scale value of this 
function by averaging it from r p — 5 to 20/i~ 1 Mpc. This 
is what we define as the relative bias of the two 

samples. 

We compute e rrors by jackknife r esampling of the data 
on the sky fe.g.. IZehavi et afll2005|) . We divide our sam- 
ple area into twenty, roughly equal-area (~ 175°), con- 
tiguous regions, and then make our measurements twenty 
times: each time dropping a different region from our 
samples. Since we do not attempt to fit a model to our 
measurements, it is not necessary to compute the full 
covariance matrix. We only calculate the diagonal er- 
rors and show them in Figures ^ and [3 These errors are 

given by a y = Eili(y* ~ v) 2 > where N = 20, and 

Di are a set of w p (r p ) or measurements made in 

the twenty jackknife samples. 

3. THE CLUSTERING OF GROUPS AS A FUNCTION OF 

MASS 

Figure Q shows the group-galaxy projected cross- 
correlation functions for our four mass-threshold group 
samples. On small scales, this essentially counts pairs be- 
tween each group center and its member galaxies. The 
cross-correlation function on these scales thus simply re- 
flects the average radial profile of galaxies within groups. 
On scales much larger than the typical size of groups, the 
cross-correlation function counts pairs between groups 
and galaxies in other groups. The transition between 
these "1-group" and "2-group" terms (analogous to the 
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Fig. 1. — Projected group-galaxy cross-correlation functions for 
groups of different mass thresholds. Gr oup masses are estimated 
from their abundances as described in § 12.21 The mass thresholds 
are listed in the panel (in units of /i _1 Mq) and the vertical arrows 
denote their virial radii. The galaxies used in the cross-correlation 
have absolute magnitudes Mo.i r < —20. Errorbars are computed 
by jackknife resampling of the data set. 



1- and 2-halo terms in the halo model) is striking in Fig- 
ure ^ much more so than in the galaxy-ga laxy autocor- 
relatio n function, where it was detected bv lZehavi et al.l 
(2004|). 

Comparing the different group mass samples on small 
scales, we see that they have cross-correlation functions 
of very similar shape, but more massive groups have a 
higher amplitude and are shifted to larger scales than less 
massive groups. Groups of different masses thus have 
radial profiles with similar functional forms, but more 
massive groups contain a higher overall density of galax- 
ies and they extend to larger radii, as expected. The four 
vertical arrows in Figure^show the virial radii of our four 
mass thresholds, estimated as R v i r = (3M/8007TP) 1 / 3 , 
where M are the thresholds and p is the mean density 
of the universe (assuming Cl m ~ 0.3). The 1- to 2-group 
transitions occur at roughly the virial radii of the mass 
thresholds, confirming that (1) the location of the tran- 
sition reflects the typical size of groups being considered, 
and (2) our estimated group masses yield virial radii that 
are in agreement with the physical sizes of the groups. 

On very small scales (< 0.3/i _1 Mpc) the cross- 
correlation functions flatten, suggesting that the radial 
profiles of groups have a core. However, on these scales 
the cross-correlation functions are very sensitive to the 
definition of group center. Any departure of our esti- 
mated centers from the "true" group centers would result 
in just this sort of flattening. 

On large scales, the relative amplitudes of the cross- 
correlation functions for the different mass groups show 
that more massive groups are more clustered (i.e., have 
a higher bias) than less massive groups, as expected. We 
quantify this dependence in Figure [21 where we show 
the bias of each sample as a function of its mass thresh- 
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Fig. 2. — Large-scale bias of groups as a function of group mass 
threshold, normalized to the bias of groups with mass greater than 
10 13 /i — 1 Mq. Data points are computed by taking the projected 
galaxy-group cross-correlation function of each group mass sam- 
ple, dividing by that of the \og(M/h~ 1 MQ) > 13 sample, and 
averaging from 5 to 20h~ x Mpc. Errorbars show the 1 — a jack- 
knife uncertainties. The two sets of points show results for group 
masses estimated using two different cosmological models, as listed 
in the panel. The two curves show the large-scale bias of dark mat- 
ter h alos as a function of mass (as computed by Scliak & Warren 
2004) for the same cosmological models. 



old (filled black points). We normalize the bias to that 
for the \og(M/h~ 1 Mo) > 13 sample and thus compute 
the bias ratio b(M)/b(M = I 13 /i _1 M Q ), as described in 
S 12.31 We obtain errorbars by jackknife resampling of the 
data on the sky. For com parison. FigureElalso shows the 
iSeliak &: WarrerJ l)2004|1 theoretical bias function of dark 
matter halos, which we also normalize to M = 13h~ 1 M Q 
(solid black curve). To first order, our measurements 
agree with the theory. The large-scale clustering ampli- 
tude of groups rises with mass, following the well-known 
halo bias function llCole fc Kaiserl Il989t IMo fc White! 
IToMlSheth fc Tormenlll99a ISeliak & Warrenll200l ~ 
In detail, however, Figure |21 shows that the bias func- 
tion for our groups rises slightly more steeply with mass 
than the halo bias function b(M). ISeliak &: WarrerJ 
(2004) showed that when halo mass is scaled by the 
nonlinear mass M*, the bias function b(M/M*) is fairly 
invariant to the values of cosmological parameters, at 
least for the limited range of parameter space explored 
by those authors (deviations from this were shown by 
iTinker et al J 120051. In other words, most of the effect 
of changing the cosmological model comes from the cor- 
responding change in M*. This means that, in order to 
have a steeper halo bias function, we must lower M*. 
Lowering M* will make two fixed masses (e.g., 10 13 and 
10 14 /i -1 M©) move to higher values of M/M*, which will 
yield a steeper bias ratio between the two masses. The 
cosmological parameter that has the strongest impact on 
is the power spectrum amplitude a%, with the mass 
density O m having a secondary effect and the spectral 
index n s placing third. The red dotted curve in Figure |21 
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shows the resulting halo bias function when as is low- 
ered from 0.9 to 0.8. As expected, the function becomes 
steeper. In order to be self-consistent, we must now also 
use this lower value of as when we assign group masses. 
Lowering as decreases the abundance of massive halos, 
which leads to lower masses for all groups. This, in turn, 
will make our fixed-mass threshold samples contain fewer 
and more rare groups (i.e., higher M/M*) that are more 
strongly clustered. We calculate new group masses using 
a s = 0.8 and re-make our samples. We show the re- 
sulting bias function in Figure [5] (open red points). The 
observed group bias function becomes slightly steeper, 
but the change is well within the errorbars and smaller 
than the corresponding change in the theoretical halo 
bias curve, especially for our lowest mass sample. The 
theory now gives a better match to the data. 

These results demonstrate that, in principle, the com- 
bination of group abundances and clustering (one to get 
the mass scale and the other to compare to theory) has 
the power to con strain cosmological parameters (see, e.g., 
IMo et al]ll99fih . For example, if we were to take our 
errors at face value and assume no systematic errors, 
our data would rule out the (Q m , erg, n s ) = (0.3, 0.9, 1.0) 
model at very high confidence. Moreover, lowering Q m 
to 0.25 or n s to 0.95 would move in the right direction, 
but would not be sufficient to give good agreement. a s 
would have to be lowered as well. However, this exercise 
is not very meaningful without carefully taking into ac- 
count the main systematic effects present in the analysis: 

(1) our groups do not perfectly correspond to halos, and 

(2) we are ignoring the scatter in the luminosity-mass 
relation when we assign group masses. These issues are 
especially important given that almost all of the con- 
straining power that we have just discussed comes from 
our lowest mass sample, where most groups only contain 
a couple galaxies. Properly addressing our systematic 
errors would require analysis of realistic mock catalogs 
constructed using a range of input cosmological param- 
eters, in order to understand the biases resulting from 
group identification, and it is beyond the scope of this 
paper. The main thing we wish to take away from Fig- 
ure [5] is that the bias of our groups agrees with that 
expected for halos. This acts as a good sanity check that 
our group samples are not too different from their under- 
lying halos. Furthermore, our data seem to prefer a low 
value for M* , in agre ement with current CMB constraints 
ijSpergel et al.ll200fil) . but we cannot say more than this 
without a detailed treatment of systematic effects. 

4. THE CLUSTERING OF GROUPS AS A FUNCTION OF 
OTHER PROPERTIES 

We have established that the large-scale clustering of 
our groups depends on mass in the manner expected of 
dark matter halos. We now move on to explore the de- 
pendence of clustering on other group properties. Since 
many group properties are likely to correlate with mass, 
it is possible that any clustering dependence we find is 
simply due to this correlation. Therefore, it is important 
that we control for mass and look for residual dependen- 
cies on other group properties at fixed mass. For each 
mass threshold sample and for each group property, we 
split the groups into a "high" half and a "low" half ac- 
cording to that property, in a way that keeps the mass 
distributions of the two halves equal. 



The properties we consider are: 

(a) (g — r)tot (Total group color): We compute a total 
color for each group by adding up all the r-band light to 
get tot, adding up all the g-band light to get M ffito t, 
and setting (g— r)tot = M S) tot — M r ,tot- This is essentially 
a luminosity-weighted color. Redder groups make it into 
the "high" samples and bluer groups make it into the 
"low" samples. 

(b) (g — r) cen (Central galaxy color): We define the 
central galaxy to be the most luminous in each group 
(in the r-band) and use its °'\g — r) color. Groups with 
redder central galaxies make it into the "high" samples 
and groups with bluer central galaxies make it into the 
"low" samples. 

(c) (g — r) avg (Average galaxy color): We take the av- 
erage °'\g — r) color of all the galaxies in each group. 
This weights all galaxies equally and thus counts satellite 
(non-central) galaxies more than the total group color. 
Again, redder groups make it into the "high" samples 
and bluer groups make it into the "low" samples. 

(d) M r ,cen (Central galaxy luminosity): This is just the 
absolute r-band magnitude of the most luminous galaxy 
in each group. In cases where iV = 1, M riCen is naturally 
equal to M r ,tot- Groups with a more luminous central 
galaxy make it into the "high" samples and groups with 
a less luminous central galaxy make it into the "low" 
samples. 

(e) N (Group multiplicity): The number of galaxies 
with Mo-v < —19.9 in each group. We consider all val- 
ues, from N = 1 isolated galaxies, to high N rich clus- 
ters. Groups with more members make it into the "high" 
samples and groups with fewer members make it into the 
"low" samples. 

(f) n ccn (Central galaxy concentration): A seeing- 
convolv ed Sersic model is fit to the i-band radial light 
profile l|Blanton et alJl2003b(l . The Sersic index n is a 
parameter of the fit and it is a measure of the concentra- 
tion of the light, n = 1 and 4 correspond to exponential 
and de Vaucouleurs profiles, respectively. Values greater 
than n ~ 5.9 are capped at that value. We use the Sersic 
index of the most luminous galaxy in each group. Groups 
with a more concentrated central galaxy make it into the 
"high" samples and groups with a less concentrated cen- 
tral galaxy make it into the "low" samples. 

(s) a v (Group velocity dispersion): This is simply 
the redshift dispers ion of each group, as described in 
iBerlind et all <|200fl) . Groups with N = 1 have a v = 0. 
Groups with a higher velocity dispersion make it into the 
"high" samples and groups with a lower velocity disper- 
sion make it into the "low" samples. 

(h) M r ,2— l (Luminosity gap): This is the magnitude 
difference between the first and second most luminous 
galaxies in each group (in the r-band): M r ,2-i = M r ,2 — 
Mr )Ce n- In groups with only one galaxy, we use the lower 
limit for the luminosity gap, which is the magnitude dif- 
ference between that galaxy and the magnitude limit of 
our sample: M r ^-i(N = 1) = 19.9 — M r ,cen- Groups 
with a larger gap make it into the "high" samples and 
groups with a s maller gap make it in to the "low" samples. 
As described in lBerlind et al.l l|2006f) . galaxies that do not 
have measured redshifts due to fiber collisions were in- 
cluded in the galaxy sample and assigned the magnitude 
and redshift of their nearest neighbor. As a result, in 
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some groups, the first and second most luminous galaxies 
have the same luminosity. Since we have no information 
about the luminosity gap in these cases, we randomly 
place half of them in the "high" samples and the other 
half in the "low" samples. 

As mentioned above, we wish to study the dependence 
of large-scale clustering on these properties at fixed mass 
in order to remove any clustering dependence coming 
from their correlation with mass. We do this as follows: 
for each candidate group, we define a mass bin of width 
AlogAf = 0.2 (our results are not sensitive to the choice 
of bin width) centered on the mass of the group, and 
we create a list of all groups whose estimated masses 
lie within the bin. We then rank this list according to 
a given group property and ask whether our candidate 
group sits in the top or bottom 50% of the ranked list. If 
it is in the top 50% we place it into the "high" sample, 
and if it is in the bottom 50% we place it into the "low" 
sample. In this way, we create two samples of roughly 
equal size that have identical mass (total luminosity) dis- 
tributions, but quite disjoint distributions in the group 
property used to split them. We do this separately for 
each of the properties in the above list. 

Figure [3] shows the result of this splitting procedure. 
Each panel shows the relation between a group prop- 
erty and M r tot (our proxy for mass). Groups that were 
placed in the "high" and "low" samples are represented 
by red and blue dots, respectively. The four vertical dot- 
ted lines show the group luminosities that correspond to 
our four mass thresholds. The boundary between the 
two samples in each case shows the median relation of 
each group property with total group luminosity. Some 
properties, like central galaxy luminosity, group multi- 
plicity, and velocity dispersion, have strong correlations 
with group luminosity for all M r ^ ot , whereas others, such 
as the three colors, central concentration, and luminosity 
gap, correlate with group luminosity only at low M r ^ot- 
The three properties that correlate the most with M r ^ t 
are simple to understand. M r ^ ccn correlates with M r ^ ot 
because the total group luminosity includes the central 
galaxy luminosity and is often dominated by it in small 
groups. At fixed M rjto t, however, we have no reason 
to think that M r ^ ccn will correlate with mass. In other 
words, M r ^ot and M riCen probably do not correlate with 
mass independently of each other. The same is not true 
for group multiplicity or velocity dispersion. These prop- 
erties likely correlate with group mass more indepen- 
dently of M r .tot- For example, at fixed group luminos- 
ity, groups with higher velocity dispersion probably have 
larger masses on average than groups with lower veloc- 
ity dispersion. Another way to say this is that a v and 
N probably correlate with the scatter in mass at fixed 

Mr, tot- 

Now that we have split the groups into halves accord- 
ing to each group property, we can measure the bias ra- 
tio 6iow/&high ( as described in § I2.3|l for each property 
as a function of group mass. The results are shown in 
Figure 0] Each type of point corresponds to a particular 
group property, as listed in the panel. Values of &i ow /^high 
equal to unity indicate that the large-scale clustering of 
groups does not depend on that group property. Values 
above unity indicate that groups in the "low" sample 
cluster more strongly than those in the "high" sample, 
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Fig. 3. — Group properties as a function of total absolute mag- 
nitude, showing the divisions into two samples for each property. 
Each panel shows a particular group property: (a) total group 
g — r color, (b) brightest galaxy g — r color, (c) mean galaxy g — r 
color, (d) brightest galaxy r-band absolute magnitude, (e) group 
multiplicity, (f) brightest galaxy sersic index, (g) group velocity 
dispersion, (h) absolute magnitude gap between first and second 
brightest galaxy. For each property, groups are placed into a "top 
50%" (red dots) or "bottom 50%" (blue dots) sample depending 
on whether their value is higher or lower than the median value 
for all groups at fixed group mass. Vertical dotted lines show the 
total absolute magnitudes corresponding to our four group mass 
thresholds. 



and the opposite is true for values below unity. 

In order to assess the significance of our results, we 
create 200 realizations of random "high" and "low" sam- 
ples, where the group property used in the splitting is 
a random number. The dark and light shaded regions 
in Figure 01 enclose 68% and 95% (i.e., la and 2a) of 
these 200 realizations, respectively and are centered on 
the median. The shaded band grows with mass because 
the higher mass samples contain fewer groups and their 
clustering measurements are therefore noisier. We have 
also computed jackknife errors for our &i ow /&high mea- 
surements and verified that they are roughly in agree- 
ment with the random splitting errors. 
We draw several conclusions from Figure 01 
(1) Only group multiplicity and velocity dispersion cor- 
relate significantly with large-scale clustering at all esti- 
mated group masses. Furthermore, they behave in ex- 
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Fig. 4. — Ratio of large-scale bias of groups with low values 
of a given property over that of groups with high values of that 
property, as a function of group mass threshold. At each mass 
threshold, different points show this bias ratio for groups split by 
different properties. Point types and properties are listed in the 
panel. All points are computed by taking the projected galaxy- 
group cross-correlation function of the "low" sample, dividing by 
that of the "high" sample, and averaging from 5 to 20h _1 Mpc. The 
shaded regions enclose 68% (dark shaded region) and 95% (light 
shaded region) of the values of £>iow/^high that result from randomly 
splitting the group sample into two halves 200 independent times, 
and thus represent 1— and 2 — a detection levels. 



actly the way that is expected if these group properties 
correlate with mass at fixed group luminosity. Groups 
with higher N or a v cluster more strongly than groups 
with lower N or a v . Moreover, the bias ratio &iow/frhigh 
is approximately the same (~ 0.7) for both properties 
and for all masses. If N and a v correlate with mass at 
fixed group luminosity, as we expect, then their bias ratio 
will depend on two factors: the amount of scatter in the 
luminosity-mass relation (after all, if the scatter is zero 
then TV and a v cannot provide additional information on 
the mass), and the steepness of the bias function b(M) 
(the steeper the function, the more a given amount of 
scatter can affect the bias) . These two effects have oppo- 
site dependencies on mass: as mass increases, the scatter 
in the luminosity-mass relation should decrease, while 
the bias function steepens. Figure 0] suggests that these 
two effects balance each other out to produce a roughly 
constant &iow/&high ratio as a function of mass. 

(2) The three group colors that we consider show evi- 
dence of a correlation with large-scale clustering at high 
mass, but not at low mass. Of these, the most signifi- 
cant is the central galaxy color. High-mass groups and 
clusters with blue central galaxies are more strongly bi- 
ased than those with red central galaxies. Actually, at 
these masses almost all central galaxies are red (see the 
second panel in Fig. so it is more accurate to say 
that groups with less red central galaxies are more clus- 
tered than those with redder central galaxies. This is a 
~ 2a effect for the 10 13 and 1O 13 - 5 /i _¥ M mass thresh- 
olds and grows to a ~ 3a effect at 10 14 /i _1 Af Q . There 



is a similar effect for total and average group colors in 
the sense that less red massive groups are more biased 
than redder groups, but it is only at the ~ 1.5a level of 
significance. Moreover, group color correlates with cen- 
tral galaxy color by construction, so it is not clear how 
much of the signal seen for total and average group color 
actually comes from the color of the central galaxy. 

(3) There is a hint that, at high masses, groups with 
more concentrated central galaxies are more biased than 
groups with less concentrated central galaxies. This may 
seem to be in conflict with the result for central galaxy 
color because usually more concentrated galaxies are red- 
der. However, at these very high luminosities, no such 
correlation exists. In any case, this is only a la measure- 
ment, so we do not consider it significant. 

(4) The luminosity gap statistic shows that groups with 
smaller gaps are more biased at low masses and less bi- 
ased at high masses. At high masses this is only a la ef- 
fect, but at the lowest mass threshold it is a ~ 3a effect. 
We must be cautious, however, because at these lowest 
masses, a large fraction of groups have N = 1, where the 
gap statistic is highly unreliable. There could therefore 
be a large systematic error that affects this point. We 
test this by repeating the measurement after throwing 
out all N = 1 groups, and we find that the detection 
disappears completely. The overall significance of the 
luminosity gap detections are therefore low. 

(5) Central galaxy luminosity does not appreciably cor- 
relate with large-scale clustering at fixed total group lu- 
minosity. Since we are comparing samples at fixed total 
group luminosity, splitting by central galaxy luminosity 
is the same as splitting by the ratio of central-to-total 
luminosity. This ratio is similar to the luminosity gap 
statistic in the sense that groups with a high central- 
to-total luminosity ratio are likely to also have a large 
luminosity gap. However, the central-to-total ratio is 
more robust than the luminosity gap because it is less 
affected by systematic effects due to incompleteness and 
group misidentification. The different behavior exhibited 
by these two statistics provides additional evidence that 
the luminosity gap results are not significant. 

Aside from the observed clustering dependence on 
group multiplicity and velocity dispersion, which are 
most likely due to their correlation with mass, our most 
significant detection is that group color (especially the 
central galaxy color) correlates with large-scale bias in 
high luminosity groups. The simplest explanation for 
this effect is that group color correlates with mass at 
fixed group luminosity, much like multiplicity and veloc- 
ity dispersion. One potential problem with this expla- 
nation is that the observed bias ratio for central galaxy 
color is more extreme than that for velocity dispersion, 
which implies that color must correlate more strongly 
with mass than does a v . This definitely seems odd, but 
perhaps it is true. Since we know that velocity dispersion 
should correlate with mass at fixed group luminosity, we 
can check whether groups with less red central galaxies 
have a different a v distribution than groups with redder 
central galaxies. 

Figure [5] shows the velocity dispersion distributions for 
the "high" (solid red histograms) and "low" (dotted blue 
histograms) samples for each group property in the case 
of log(M I 7i _1 Mq) > 14 groups. Each panel also lists 
the mean values of a v for the two samples. With the ex- 
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Fig. 5. — Velocity dispersion histograms for samples of mas- 
sive (M > 10 14 h~ 1 Mq) groups split by different group properties. 
Each panel shows a particular group property, as in Fig. [3] In each 
panel, the solid and dotted histograms show the velocity dispersion 
distribution for groups with high and low values of the property, re- 
spectively. The mean velocity dispersion for the "high" and "low" 
samples are denoted {<j v )a and {ct v )b, respectively, and are listed 
in each panel. 



ception of multiplicity and, of course, velocity dispersion 
itself, all other properties have "high" and "low" sam- 
ples with nearly identical distributions. In other words, 
these properties show no correlation with o~ v . This result 
provides strong evidence that whatever trends exist in 
Figure 01 are not due to correlations with mass. The one 
notable exception is group multiplicity, which correlates 
with a v , as expected. 
We discuss the implications of these results in § 

5. THE CLUSTERING OF HALOS 

In the previous section, we showed that there is a sig- 
nificant dependence of high mass group clustering on a 
group property other than mass, and that this depen- 
dence is most likely not due to a correlation with mass. 
One possibility is that this is a detection of the recently 
discovered assembly bias for dark matter halos; namely, 
that halo clustering depends on halo age/concentration 



at fixed mass. iWechsler etTafl (|2006|) found that for 
M > M*, low concentration halos are more strongly bi- 
ased than high concentration halos. However, the simu- 
lations used in that study were fairly small in volume (al- 
beit very high in resolution) making it difficult to study 
the high mass regime whe re halos are rare. By necessity, 
the IWechsler et all (|2006j) results for high M/M* came 
from high redshift outputs of the simulation, where M* 
is lower. IWetzel et al.l l|2006f) investiga ted this with a 
much larger simulation and confirmed the We chsler et alJ 
(2006) results for high mass halos. 

We wish to compare the age/concentration effect for 
high mass halos to our results for galaxy groups. We 
therefore need to know the ratio 6i ow /&high for halos, 
where halos are split into "low" and "high" samples in a 
similar way as groups were split in § We do this using 
a set of N-body simulations, which we describe below. 

5.1. N-body Simulations and Analysis 

We use 20 independent N-body simulations of a ACDM 
cosmological model, with O m = 0.27, ft a — 0.73, 
n b = 0.046, h = flo/(100 kms" 1 Mpc" 1 ) = 0.72, 
n s = 1.0, and erg = 0.9. Each simulation follows 
the evolution of 512 3 dark matter particles of mass 
m p = 7.49 x lO lo /i _1 Af0 in a periodic box of length 
512/i _1 Mpc. Initial conditions were set up using second- 
order Lag rangian perturbatio n theory (2LPT) as de- 
scribed bv lCrocce et alJ l|2006D , an d the simulatio ns were 
run using the GADGET2 code ijSpringell 12005^ . The 
gravitational force softening was set to e grav = 40/i _1 kpc 
(Plummer equ ivalent). The simula tions are described in 
more detail in lCrocce et ]d](|2006). 

We identify halos in the z — dark matter particle 
distributions using a friends-of-friends algorithm with a 
linking length equal to 0.2 times the mean inter-particle 
separation. The mass of each halo is set to the sum 
of its me mber particle masse s, times a correction factor 
given bv lWarren et alJ l)2006l) that accounts for system- 
atic effects due to the finite number of particles per halo. 
For the large halo masses that we consider in this paper, 
this correction is never larger than 1.5%. This procedure 
yields a total of 90231 halos (in all 20 simulations) of 
mass greater than 10 14 ft. -1 Mq. For each halo, we com- 

/ \ 1/3 

pute a virial radius equal to i? V j r = ( sooxp ) ' wnere 
M is the halo mass and p is the mean density of the uni- 
verse. Finally, we assign center-of-mass coordinates to 
all halos. 

Halo concentrations are simpler to measure than ages 
because they can be measured from the z = simula- 
tion outputs alone and do not require the construction 
of merger trees. For this reason, we use concentrations in 
this study; however, we note that similar results should 
hold for formation ti me since these halo pr operties corre- 
late with each other ( Wechs ler et al.ll2"f)02D . We measure 
a concentration for each halo in the simplest possible 
way: we choose a fixed fraction of the halo virial radius 
and measure the ratio of mass enclosed within this ra- 
dius to the total mass. We choose a value for this fraction 
equal to 0.3, but we also compute concentrations using 
values as low as 0.1 and as high as 0.5 in order to test 
how robust our results are to this definition. 

We then split all halos of mass greater than 
10 14 /i _1 Mq into "high" and "low" samples, in exactly 
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the same way as we split our groups in § 0] For each 
halo, we create a mass bin of width AlogM = 0.1 that 
is centered on that halo's mass, and we make a list of 
all halos that lie within this bin. We then sort this list 
by concentration, and we place the halo in question in 
the "high" or "low" sample depending on whether its 
concentration places it in the top or bottom 50% of the 
sorted list. We also create a second set of samples by 
choosing halos whose concentrations place them in the 
top or bottom 25% of halos within their mass bin. In 
this way, we create halo samples of high and low con- 
centration that have the same mass distributions, thus 
controllin g for the well-k nown concentration-mass rela- 
tion jBullock et alJl200l|) . 

We measure halo-halo autocorrelation functions for the 
"high" and "low" samples in each of the 20 simulations. 
We use the simple estimator £(r) = DD/RR — 1, where 
DD is the number of halo-halo pairs, and RR is the 
number of random-random pairs, which we calculate an- 
alytically as RR = (l/2)(4/3)7rn 2 y b ox(r3 ut - rfj, where 
n is the mean density of halos, Vbox is the volume of the 
simulation cube, and r ou t and r; n are the outer and inner 
radii of the bin in which £(r) is being calculated. We then 

calculate bias functions b{r) = ^/6owM/6iighM- We 
compute average correlation functions and bias functions 
from the 20 simulations, and we calculate errors from the 
dispersion among these independent realizations. 

5.2. Results 

The top panel of Figure HJI shows halo correlation func- 
tions for the 50% "low" and "high" samples (thick red 
curves), as well as the 25% "low" and "high" samples 
(thin blue curves) . All correlation functions show the ex- 
pected turnover at low scales due to the fact that friends- 
of- friends halos cannot overlap. There can thus be no 
halo pairs at separations less than twice the virial ra- 
dius of the smallest halos considered. The bottom panel 
shows low/high bias functions for the 50% (thick red 
curve) and 25% (thin blue curve) samples. 

Figure clearly shows that low concentration halos 
have a significantly higher amplitude of £(r) than high 
concentration halos, at fixed mas s. This is a very high 
sign al-to-noise confirma tion of the lWechsler eTall ( 2006) 
and lWetzel et alJ |2006j) results. The bias functions be- 
come scale- independent at large scales, showing that halo 
concentration affects the amplitude, but not the shape of 
the correlation function, as expected. The low/high bias 
of the 25% samples is naturally higher than that of the 
50% samples, but it is interesting that all the difference 
comes from the high concentration end. The correlation 
functions of the lowest 50% and 25% concentrated halos 
have exactly the same amplitude. This means that the 
relation between concentration and halo bias, at fixed 
mass, is flat below the median concentration, and only 
starts dropping at higher concentrations. 

We compute large-scale asymptotic values of 6i ow /^high 
by averaging the bias functions in the range where they 
are scale independent: from 9 to 35/i -1 Mpc. We do 
this separately for each of the 20 simulations and then 
compute the mean and uncertainty in the mean, which 
we calculate from the dispersion among the simulations. 
The relative bias of low over high concentration halos 
is friow/frhigh = 1-243 ± 0.017 for the 50% samples and 
&iow/frhigh = 1.316 ± 0.032 for the 25% samples. These 
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Fig. 6. — Top panel: Two-point autocorrelation functions for 
simulated dark matter halos of mass greater than 10 14 /i — 1 Mq, 
split into subsamples by concentration. Thick red solid and dashed 
curves show correlation functions of the 50% least and most concen- 
trated halos, respectively. Thin blue solid and dashed curves show 
the same for the 25% least and most concentrated halos. Bottom 
panel: Bias function defined as the square root of the correlation 
function for low concentration halos divided by that for high con- 
centration halos. The thick red, and thin blue curves shows the 
bias for the 50% and 25% concentration splits, respectively. Hori- 
zontal doted lines show the values of the bias, averaged from 9 to 
35/i Mpc in these two cases. The values of this large-scale bias 
are 1.243 ± 0.017 and 1.316 ± 0.032 for the 50% and 25% concen- 
tration splits, respectively. In both panels, results are averaged 
over 20 independent N-body simulations, and errors show the un- 
certainty in the mean estimated from the dispersion among these 
simulations. 



values are shown in Figure El as dotted lines. We note 
that our results are not sensitive to the specific defini- 
tion of concentration that we use. When we measure 
concentrations using different fractions of the virial ra- 
dius, the values of &k>w/&high change by less than their 
uncertainties. In the next section, we discuss the connec- 
tion between these results for halos and those for groups 
shown in the previous section. 

6. DISCUSSION 

In § 0] we found that the large-scale bias of massive 
groups correlates with their central galaxy color, at fixed 
estimated group mass. Specifically, massive groups with 
less red central (brightest) galaxies cluster more strongly 
than those with redder central galaxies. Furthermore, 
we showed that this trend is probably not due to a cor- 
relation of central galaxy color with true mass at fixed 
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estimated mass. In other words, we have found that 
the large-scale clustering of massive groups depends on 
a "second parameter" other than mass. This result can 
be understood i n the cont ext of the recent fin dings by 
iWechsler et aD Ipffll and iWetzel eTal! IpOOl . as well 
as our results in § El showing that massive dark mat- 
ter halos also have a "second parameter" , whether it be 
concentration, age, or some other feature of the halo as- 
sembly history. If central galaxy color correlates with 
halo concentration or age at fixed halo mass, then it will 
also correlate with large-scale bias, since halo concentra- 
tion and age correlate with bias. Moreover, the effect 
should only be present at high masses, and should dis- 
appear as M approaches Af*, which is exactly what we 
find. 3 This means that our result for groups most likely 
constitutes observational evidence for the recently dis- 
covered "second parameter" effect for halos. The only 
alternative is that the color of a halo's central galaxy 
somehow "knows" about the density field at large scales 
(~ 10/i _1 Mpc), independently of its halo's assembly his- 
tory. Though possible in principle, this is highly unlikely. 

For groups of mass greater than lO 14 /i _1 M0, we found 
that the relative bias between groups with the bottom 
and top 50% of central galaxy color is &i ow /&high ~ 
1.5 ±0.2. 4 For comparison, the relative bias between the 
50% low and high concentration halos in our simulations 
is friow/frhigh = 1-243 ± 0.017. This bias ratio is closer to 
what we found for total and average group color. This 
difference could indicate that central galaxy color corre- 
lates with a feature of the halo assembly history that is 
more directly tied to large-scale bias than concentration. 
However, the difference between these numbers is not 
highly statistically significant. This highlights the need 
for repeating these observational investigations with a 
much larger sample volume that contains a much larger 
number of massive groups. 

Our results imply that massive halos that formed ear- 
lier contain redder central galaxies than halos of the same 
mass that assembled more recently. This is interesting 
because it is not obvious what one should expect. One 
could argue that recently assembled halos are more likely 
to have had a recent major merger, which might result 
in a central galaxy with no star formation, whereas older 
halos might contain central galaxies that have had time 
since their last merger to accrete gas and have some on- 
going star formation. On the other hand, one could argue 
that in this high mass regime, there will be no cool gas to 
accrete and that older halos will simply have older cen- 
tral galaxy stellar populations. Our results suggest the 
latter scenario. 

It is interesting that the luminosity gap, which we ex- 
pected to be the group property most directly related 
to age, showed no significant correlation with large-scale 
bias at fixed mass. The reasoning behind our expecta- 
tion is that in old groups, satellite galaxies will have had 
more time to merge with the central galaxy, resulting in 
a very massive central galaxy and thus a large luminos- 
ity gap. If anything, our results show the opposite trend 
than expected: massive groups with a high luminosity 
gap (and hence older) are slightly more clustered than 

3 The effect should re-appear with the opposite sign at lower 
masses, but our group catalog does not probe that mass regime. 

4 The uncertainty comes from jackknife resampling. 



groups with a low gap. This contradicts the theoretical 
results for halos. One possible explanation for this is that 
the luminosity gap does not actually correlate with age 
in the way we expect. However, it could also be that our 
measured luminosity gaps suffer severely from system- 
atic effects due to incompleteness and group misidentifi- 
cation and are, therefore, not representative of the true 
luminosity gaps within halos. We believe this latter ex- 
planation because we find quite different results when we 
study the dependence of clustering on central galaxy lu- 
minosity, which should act like the luminosity gap. Old 
groups should have both a larger luminosity gap, and a 
more luminous central galaxy than young groups of the 
s ame mass. 

Zcnt ner et al.l <|2005|) predicted that halo concentra- 
tions and formation times should be correlated with the 
number of subhalos within the halo in the sense that 
older, more concentrated halos have fewer subhalos. This 
correlation arises because older halos accrete their sub- 
halos earlier and thus give them more time to sink to 
the center and be destroyed. If this prediction is cor- 
rect, we should see a trend whereby high mass groups 
with many member galaxies (high N) are more strongly 
biased than poorer groups of the same mass. Moreover, 
this trend should disappear at smaller group masses since 
the bias-age trend also goes away at smaller masses. In 
this paper, we find the predicted trend at high masses, 
but, unlike the prediction, it persists at lower masses. It 
is n ot straightforward to interpret our result in terms of 
the iZentner et alJ l|2005fl prediction because we do not 
compare groups at fixed mass, but rather at fixed total 
luminosity. If group multiplicity N correlates with group 
mass at fixed group luminosity then we also expect a 
trend like the one we find. Since the trend we see per- 
sists at low group masses, we conclude that we are likely 
detecting a multiplicity-mass correlation rather than an 
i ntrinsic mult i plicity -bias relation. 

lYang et alJ 1)2006(1 have also detected a "second param- 
eter" effect for groups. In a study similar to ours, but 
using 2dFGRS groups, they found that groups contain- 
ing central galaxies with low SFR are more biased than 
groups containing central galaxies with high SFR. More- 
over, they found this effect at all group masses. This 
seems to be at odds with our results because if low SFR 
impli es redder color, th en we h ave found the opposite ef- 
fect. lYang et all l)2006|) use the lMadewick et alJ l|20021 r\ 
parameter as a proxy for SFR, which reflects the average 
emission- and absorption-line strength in the galaxy rest 
frame spectrum. This is not necessarily correlated with 
g — r color for the luminous red galaxies that are the 
central galaxies of these massive systems. Nevertheless, 
at its face, this result looks completely opposite to ours 
and remains an interesting puzzle. 

7. SUMMARY 

In this paper we have investigated the clustering of 
galaxy groups in the SDSS. The principal goal of this 
study was to look for secondary dependences of large- 
scale clustering on group properties other than mass. We 
have estimated group masses from their abundances us- 
ing the total group luminosity as a proxy for mass. Our 
main results are: 

1. The measured large-scale bias of groups as a func- 
tion of estimated mass is in agreement with the theoret- 
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ical halo bias function, given a standard ACDM cosmo- 
logical model. The measurements suggest a preference 
for a low value of erg, in agreement with current CMB 
constraints, but this result is not robust since it could be 
subject to systematic effects. 

2. We have measured the residual dependence of 
group bias on other group properties, at fixed estimated 
mass. The properties we have considered are: total group 
color, central galaxy color, average group color, central 
galaxy luminosity, group multiplicity (richness), central 
galaxy concentration, group velocity dispersion, luminos- 
ity gap between first and second brightest galaxies. Of 
these properties, only group multiplicity, velocity disper- 
sion, and central galaxy color show a significant corre- 
lation with bias at fixed estimated mass. Specifically, 
groups with higher multiplicity, higher velocity disper- 
sion, and less red central galaxies cluster more strongly 
than groups with the opposite properties. The effect for 
multiplicity and velocity dispersion occurs at all masses, 
whereas the effect for central galaxy color is only signifi- 
cant at high group masses (M > 10 14 /i -1 Mq). 

3. The dependence of large-scale bias on group mul- 
tiplicity and velocity dispersion can be simply explained 
if these properties correlate with true mass at fixed esti- 
mated mass. However, the dependence of bias on central 
galaxy color cannot be explained this way and is most 
likely a real effect. This is likely observational evidence 
of recent theoretical findings that halo bias depends on 
a "second parameter" other than mass, such as age or 
concentration. 

4. Our results imply a connection between halo age 
and central galaxy color for massive halos. Halos that 
assembled earlier likely contain redder central galaxies 
than recently assembled halos of the same mass. 

5. In order to compare our results to theory, we have 
quantified the dependence of halo bias on concentration 
for high mass halos, using a set of large N-body simu- 
lations. We find that low concentration halos are more 
biased on large scales than high concentration halos, at 
fixed mass, thus confirming previous results at very high 
signal-to-noise. Specifically, we find that for halos with 
M > l0 li h~ 1 M Q , the relative bias between the 50% least 
and most concentrated halos is &i ow /&high = 1.243±0.017, 
while the relative bias between the 25% least and most 
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